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Abstract 

A nonlinear small-world network model has been presented to investigate the effect of non- 
linear interaction and time delay on the dynamic properties of small-world networks. Both 
numerical simulations and analytical analysis for networks with time delay and nonlinear inter- 
action show chaotic features in the system response when nonlinear interaction is strong enough 
or the length scale is large enough. In addition, the small- world system may behave very differ- 
ently on different scales. Time-delay parameter also has a very strong effect on properties such 
as the critical length and response time of small-world networks. 
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1 Introduction 

Since the pioneer work of Watts and Strogatz [1] on small-world networks, there arise a lot of 
interesting research on the theory and application of small- world networks [2-7]. The properties 
of complicated networks such as internet servers, power grids, forest fires and disordered porous 
media are mainly determined by the way of connections between the vertices or occupied sites. 
One limiting case is the regular network with a high degree of local clustering and a large average 
distance, while the other limiting case is the random network with negligible local clustering and 
small average distance. The small- world network is a special class of networks with a high degree of 
local clustering as well as a small average distance. Such small-world phenomenon can be obtained 
by adding randomly only a small fraction of the connections, and some common networks such as 
power grids, film stars networks and neural networks behaves like small- world networks [2-9]. 

The dynamic features such as spreading and reponse of the network have also been investigated 
in recent studies [2,3] by using shortest paths in system with sparse long-range connections in the 
frame work of small-world models. A simple time-stepping rule has been used to similuate the the 
spreading of some influence such as a forest fire, an infectious disease or a particle in percolating 
media. The influence propagates from the infected site to all uninfected sites connected to it via a 
link at each time step, whenever a long-range connection or shortcut is met, the influence is newly 
activated at the other end of the shortcut so as to simulate long-range sparkling effect such as the 
infect site (e.g., a person with influenza) suddenly travels to a new place, or a portable computer 
with virus that start to connect to the network a new site. These phenomena have been sucessfully 
studied by Newman and Watts model [2] and Moukarzel [3]. Their models are linear model in the 
sense that the governing equation is linear and the response is immediate as there is no time delay 
in their models. 

However, in reality, a spark or an infection can not start a new fire spot or new infection im- 
mediately, it usually takes some time A, called ignition time or waiting time, to start a new fire 
or infection. In addition, a fraction of infected site shall recover after a futher time of T to nor- 
mality. Thus the existing models are no longer be able to predict the response in the networks or 
systems with time delay. Furthermore, the nonlinear effect such as the competition factor as in 
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the population dynamics, congestion feature such as the traffic jam in internet communication and 
road networks, and the frictional or viscous effect in the interaction of the vertices. When consid- 
ering these nonlinear effects, the resulting small-world network model is generally no longer linear. 
Therefore, a nonlinear model is yet to be formulated. 

The main aim of this paper to present a more general nonlinear model for the small-world 
networks by extending the existing Newman- Watts [2] and Moukarzel [3] models to investigate the 
effects of time-delay, site recovery and the nonlinear interaction due to competition and congestion. 
The new model will generally lead to a nonlinear difference differential equation, whose solution 
is usually very difficult to obtain if it is not impossible. Thus the numerical simulation becomes 
essential [f]. However, we will take the analytical analysis as far as possible and compare with 
the results from numerical simulations. The characteristic chaos of the network dynamics is then 
studied by reducing the governing equation into a logistic equation. The control of the chaos is also 
investigated by introducing the negative feedback with time delay to the small-world networks. 



2 Nonlinear Model for Small- World Networks 

To investigate the nonlinear effect of time delay on the properties of a small-world network, we 
now consider a randomly connected network on a d-dimensional lattice [f ,2] (with d = 1, 2, ...), and 
overlapping on the network are a number of long-range shortcuts randomly connecting some vertices, 
and the fraction of the long-range shortcuts or probability p is relative small pCl. Now assuming 
an influence or a pollutant particle spreads with a constant velocity u = 1 in all directions and a 
newly infected site in the other end of a shortcut will start but with a time delay A. Following 
the method by Newman and watts [2] and Moukarzel [3], the total influenced volume V(t) comes 
from three contributions: one is the influenced volume with Td t d ~ 1 dt where t is time and T d is a 
shape factor, the other contribution is 2pV(t — £ — A) for a hypcrsphere started at time (. These 
two components have used studied earlier [2,3] although without the time delay parameter. Now we 
add the third component due to nonlinear interaction such as friction, slow down due to congestion 
as in the case of internet network and road traffic jam and lack of other resource as lack of oxygen 
for the fire spark to start a new fire. By assuming this nonlinear effect as —/j,V 2 (t — ( — A) where 
li <C f is a measure of nonlinear interaction. By using a continuum approach to the network, then 
V(t) satisfies the following equation with time delay 

v(t) = r d f c d - x [i + C d v(t - c - A) - ^v 2 (t - c - A)}d(, (l) 

Jo 

where d = 1,2,... and T d is shape factor of a hypersphere in d-dimensions. The Newman-Watts 
length scale [2] can be conveniently defined as 
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^ = ( P kdY/ d ' (2) 

where k — const being some fixed range. By proper rescaling t 

T = t(T d r d (d-l)l) 1/d , 5 = A(T d C d (d-iy.) 1/d - (3) 
and rewriting (1) in the rescaled form 



V(t) = j\r (f-'il + C d V(C -6)- ^V 2 (C - 5)}d(, 



(4) 



^=t d + V(T-S)-rt d V 2 (T-6), (5) 



we have a time-delay equation, after differentiating the equation d times 

d n V 

which is a nonlinear delay differential equation, whose explicit solutions is not always possible. In 
addition, the nonlinear term and time delay can have strong effect on the behaviour of the dynamic 
properties of the small- world networks. 
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3 Chaos in Small- World Networks 



From the theory of dynamical systems, it is expected that the dynamic features can be shown 
more clearly by using the representation in Poincare plane [10], which usually transform a nonlinear 
differential equation into a nonlinear iterated map or logistic equation. Now we write equation (5) 
in a difference form and take dr = S to get a logistic equation. In order to focus on the main 
characteristics of the dynamics, for simplicity, we can take 5 = 1 in 1-D (or d — 1), and we then 
have 

V n+1 =i + 2V n - ^Vl (6) 
where V n +i = V(t) and V n = V(t — 1). By changing variables 



^-V±g*<y M+ A), „„.<™>(v; + .4). A- <wf-\ (-) 

we can rewrite equation (6) as 

v n+ i = A«„(l - v n ), A = (Vl + 4/i£ 2 + l), (8) 

which is a standard form of well-known logistic equation [10]. Since the parameter range of A for 
period doubling and chaos is known, we have 

The system becomes chaotic as A is bigger than A* « 3.5699 but usually below 4.0, so the chaos 
begin at 

For A less than Aq ~ 3.0, the system approach to a fixed point, that is 



For a fixed fi, when £ < £ < £*> then A < A*, the system is in a period doubling cascade. When 
£ > £*, the system is chaotic. Clearly, as \i — > 0, — > oo. The system behavior depends on the 
lcngthscale of small- world networks. The system may looks like chaotic on a large scale greater than 
the critical length scale and the same system may be well regular on the even smaller scale. So 
the system behaves differently on different scales. 

To check the analytical results, we also simulated the scenario by using the numeric method [1,2] 
for a network size N = 500, 000, p = 0.002 and k = 2 on a 1-D lattice. Different values of the 
nonlinear interaction coefficient \i are used and the related critical length when the system of 
small-world networks becomes chaotic. Figure 1 shows for different values of \x. The solid curve 
is the analytical results (10) and the points (marked with o) arc numerical simulations. The good 
agreement verifies the analysis. However, as the typical length increases, the difference between 
these two curves becomes larger because the governing equation is main for infinite size network. So 
the difference is expected due to the finite size of the network used in the simulations. 



4 Feedback and Chaos Control of Small- World Networks 

The occurrance of the chaos feature in small- world networks is due to the nonlinear interaction term 
and time delay. This chaotic feature can be controlled by adding a negative feedback term [11,12]. 
In reality, the influence such as a signal or an influence (e.g.. influenza) only last a certain period of 
time T, then some of the influenced sites recover to normality. From the small-world model equation 
(1), we see that this add an extra term 0V(t — A — T), which means that a fraction (/?) of the infected 
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sites at a much earlier time (t — A — T) shall recover at t. So that we have the modified form of 
equation (5) as 

-^=i d + V{r-5)- K d V 2 (r -S)- K d V{r - 5 - to), (12) 

where r = T(Ydi~ d {d — l)!) 1 ^. For d = 1, we can take To = jS (j=l,2,...) without losing its 
physical importance. By using the transform (7), we have the modified form of the logistic equation 

v n+1 = Xv n (l - v n ) + a(v n - V n -j), j = 1,2, ... (13) 

with 

A = (y/1 + 4^ 2 + 1), a = ^, (14) 

which is essential the Escalona and Parmananda form [12] of OGY algorithm [11] in the chaos control 
strategy. We can also write (13) as 

v n+1 = Av n (l - v n ) - av n _j, j = 1, 2, ... (15) 

with 

A = ( v / (l-^) 2 +4/iC 2 + l), a = 0£. (16) 

This last form (15) emphases the importance of time delay and the effect of negative feedback in 
controlling the chaos. 

For a fixed value of A = 3.8, we find a critical vaule of a* = 0.27 for j = 1 and a* = 0.86 for 
j = 2 to just control the chaos so that the system settles to a fixed point. For a < a*, the feedback 
is not strong enough and the chaos is not substantially suppressed. For a > a*, the strong feedback 
essentially control the chaos of the small-world networks. Figure 2 shows that the effect of recovery 
of the infected site or the delay feedback on the system behavior. The dotted points are for the 
chaotic response when there is no feedback (A = 3.8, a = 0), while the solid curve corresponds to 
the just control of the chaos by feedback (A = 3.8, a = 0.27). This clearly indicates that the proper 
feedback due to healthy recovery and time delay can control the chaotic reponse to a stable state. 



5 Discussion 

A nonlinear small- world network model has been presented here to characterise the effect of nonlinear 
interaction, time delay, and recovery on small- world networks. Numerical simulations and analytical 
analysis for networks with time delay and nonlinear interaction show that the system reponse of the 
small- world networks may become chaotic on the scale greater than the critical length scale , and 
at the same time the system may still quite regular on the smaller scale. So the small-world system 
behaves differently on different scales. Time-delay parameter 6 has a very strong effect on properties 
such as the critical length and reponse time of the networks. 

On the other hand, in order to control the possible chaotic behaviour of small-world network, 
a proper feedback or healthy recovery of the infected sites is needed to stable the system reponse. 
For a negative delay feedback, comparison of different numerical simulations suggests that a linear 
recovery rate j3 or a linear feedback can properly control the chaos if the feedback is strong enough. 
This may has important applications in the management and control of the dynamic behaviour of 
the small- world networks. This shall be the motivation of some further studies of the dynamics of 
small- world networks. 
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Figure 2 Comparison of chaos and chaos control due to delay feedback. The dotted points are for 
the chaotic response when there is no feedback (A = 3.8, a = 0), while the solid curve corresponds 
to the just control of the chaos by feedback (A = 3.8, a = 0.27). 
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